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We present evidence for the existence of a spontaneous instability towards an orbital loop-current 
phase in a multiorbital Hubbard model for the Cu02 planes in cuprates. Contrary to the previously 
proposed 6ji phase with intra-unit cell currents, the identified instability is towards a staggered 
pattern of intertwined current loops. The orbitally resolved current pattern thereby shares its 
staggered character with the proposal of d-density wave order. The current pattern will cause 
a Fermi surface reconstruction and the opening of a pseudogap. We argue that the pseudogap 
phase with time-reversal symmetry breaking currents is susceptible to further phase transitions and 
therefore offers a route to account for axial incommensurate charge order and a polar Kerr effect in 
underdoped cuprates. 


I. INTRODUCTION 

There is now considerable evidence in underdoped 
cuprate high-temperature superconductors for a cascade 
of phase transitions, starting at high temperatures with 
the pseudogap onset at T*, followed by incommensurate 
charge order (ICO) at Too < T* and superconductivity 
at Tc < Tea- In addition, broken time-reversal symmetry 
has been associated with T* and a Kerr rotation is mea¬ 
sured below a temperature TKerr with TKen- ~ 0.75T* 
over wide doping range. At present, there is no unifying 
theory that explains this intriguing sequence of transi¬ 
tions. 

Ultrasound spectroscopy suggests that T* corresponds 
to a true thermodynamic phase transitiorf^ ; this find¬ 
ing challenges the viewpoint that the pseudogap arises 
as a correlation induced phenomenon in a symmetry un¬ 
broken paramagnetic phase. Below T* spin-polarized 
neutron scattering experiments detected weak magnetic 
moment^^lf^. These moments appear to preserve the 
translational symmetry of the lattice, and led to the pro¬ 
posal of intra-unit cell loop currents (LCs)P. However, 
the so-called 6jj LC phase by itself has difficulty explain¬ 
ing the partial gapping of charge excitation^. While 
variational methods favored the existence of LC phases 
in finite cluster^^, alternative numerical ly exa ct analy¬ 
ses reported no evidence for the On phas^^SEI] 

The Kerr effeclP^ESl that sets in below Txerr is further 
evidence for time-reversal symmetry breaking, but also 
requires that mirror symmetries be brokerP^. Through¬ 
out much of the cuprate phase diagram, TKerr and 
gtQgg^ which has motiv ated further propos¬ 
als in whic h fluctuating charge-li^l^H qj. pair-density wave 
state^^MH generate spontaneous current patterns with 
broken mirror symmetries. These scenarios assume a 
heirarchy of transitions associated with distinct symme¬ 
try breakings needed to form the fully ordered density- 



FIG. 1: (Color online) Staggered pattern of spontaneous loop 
currents. Open circle, “x” and “y” symbols denote Cud„, 2 _y 2 , 
Opx and Opy orbitals, respectively. Currents along p-d bonds 
(black arrows) are about three times stronger than those in 
p-p bonds (green arrows). 

wave state. Other proposals follow the common theme 
that the pseudogap results from the competition between 
two or more order parameter^^ZHMI 

The ICO phase involves predominantly a ch arge trans- 
fer between oxygen orbitals in the Cu 02 planes^fin^MIIll] 
This challenges notions of immutable Cu02 bands, and 
points to the necessity to employ multiorbital models for 
the ICO phastPS*^. Here, we support this reasoning and 
show that orbital resolved intra-unit cell physics is im¬ 
portant throughout the pseudogap regime. 

In this work we report the results of an unbiased cal¬ 
culation for a three-band model of Cu 02 planes which 
verifies the existence of an instability towards a stag¬ 
gered pattern of intertwined LCs (Fig. [^. This “ttLC” 
phase is different from the anticipated On phase, but 
shares its ordering wavevector Q = (tt, tt) with the ear¬ 
lier phenomenologi cal p roposal of LCs in the d-density 
wave (DDW) stat^^MIl i^j^g t^LC phase the Fermi 
surface reconstructs to form hole pockets with a concomi¬ 
tant pseudogap-like structure in the electronic spectrum. 
For realistic parameters, the ICO reported previously for 














2 


the same model HamiltoniarPH is subleading to the ttLC 
instability. Yet, the presence of staggered order favors a 
subsequent instability towards ICO with axial wavevec- 
tors connecting the tips of the hole pocket^^, consistent 
with experiments. The charge modulation of the latter 
necessarily breaks mirror symmetries and will hence al¬ 
low for a polar Kerr signaP^. This scenario is offered 
as a proposal for the cascade of phase transitions in the 
pseudogap regime of underdoped cuprates. 


II. HAMILTONIAN 

The unit cell of a single Cu02 plane is shown in Fig. 
along with the choice of orbital phases and the corre¬ 
sponding signs of the hopping terms. 

The non-interacting part of the three band model is 
given by 

Hq = Cianictcr + ^iaiP^iaa^jPcr ( 1 ) 

iota iotjlBcr 

where i and j are unit cell labels, a and /3 are orbital 
labels, (T is the spin label, eia is the orbital energy, Uiacr 
is the number operator, ti^jp is the tunneling matrix el¬ 
ement between orbital ia and j/3, and Ciaa and are 

annihilation and creation operators. Below, we suppress 
the spin labels. 

Using the translational invariance, iJp can be Fourier 
transformed to reciprocal space: 

i?o=^vi,[Ho(k)vI-k (2) 

k 

where 'E']^ = cl^, cly], and (cka) is the creation 
(annihilation) operator for an electron with crystal mo¬ 
mentum k and orbital a. Explicitly, 

Ck„ = (3) 

^ i 

T, = ( 4 ) 

^ i 

where N is the number of unit cells in the system, and 
is the position vector of a’th orbital in Fth unit cell. 
i/o(k) is readily obtained by plugging Eqs. (§ and @> 
into Eq. Q: 

( Cd ‘^iipd^y \ 

‘litpdSx '^tppSxSy I (5) 

‘2'ltpd 6y j 

where Sx = sin(fca;/2) and Sy = sin(fcy/2). A more con¬ 
venient form of Hq is obtained after the following gauge 
transformation: 


^\iX 



(6) 

cL 

— > 

-*cL 

(7) 

Cky 

-> 


(8) 

^ky 



(9) 



FIG. 2: (Color online) Unit cell of a Cu02 plane (dashed 
box). The orbital phase convention is indicated by the sign of 
the hopping matrix elements. Numbers in green enumerate 
the inequivalent bonds, and the direction of the blue arrows 
indicates a positive sign of the current flow for the current 
operator definitions given in Appendix [X] 


Hence, the final form of Hq is obtained: 


/ td ‘^tpd^x ‘^^pd^y 

Ho(k) = I ‘^^pd^x Cp AtppSxSy 

‘^tpdSy 4^tppSxSy Cp 


( 10 ) 


We set tpd = 1 so that it defines the unit of energy. 

The interacting part of the Hamiltonian includes the 
intra- (C/q) and inter-orbital (Uaj/s) Coulomb interac¬ 
tions 


iJ' = 


E 

iOL(7,j0<7' 






2Q:,j/3 




( 11 ) 

where i,j are unit cell indices, a, /3 are orbital labels and 
Viajp is nonzero for nearest-neighbors only. Throughout, 
we suppress the spin index cr, set tpp = —0.5, Cd — Cp = 
2.5, Ud = 9, Up = 3, Vpd = 2.2, and Vpp = 1 unless 
otherwise stated. The Hamiltonian H — Hq -|- H' is thus 
the conventional three-band model of cuprate^^ with a 
typical parameter selP^. 


III. INTERACTING CURRENT 
SUSCEPTIBILITY 

The current operator associated with the bond be¬ 
tween sites ia and is Jiajp = -*Uaj73(cLs73“^i/3^i«) 
where Uajp is the corresponding hopping matrix ele¬ 
ment. If (Jiajjs) > 0 then current flows from j/3 to ia. 
In momentum space, the current operator along bond m 
is given by 

im(q) = -iE[^“/3(’^’‘l)^ka4+q/3-^™a(k, q)Ck/34+qa] 

k 

( 12 ) 
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where a, /3 are the orbitals associated with the bond and 
the matrix elements of the current operators q) 

are listed in Table [T] As shown in Fig. there are eight 
distinct bonds on the Cu02 lattice. Accordingly the cur¬ 
rent susceptibility, 


= / dre“'^'^(J„^(q,r)J„(-q,0)), (13) 

do 

is an 8x8 matrix, where uJi = 2tt£T denotes the bosonic 
Matsubara frequencies. Each matrix element can be de¬ 
composed as 

Xmn Xoi^oi'Xct00'ot' ' XjBoijB'a.' 

where 

^ (k, q)/.!)^, (k',-q) 

kk' 

X f dTe*‘^'^(cj^g(T)Ck+qe-(T) 

^0 

xcL^(0)cu._qy(0)). (15) 

Previously, we investigated charge instabilities in the 
same three-band model with non-local interactions using 
a generalized random phase approximation (gRPAJpSSl 
While methods like QMC or cluster DMFT are at first 
glance more desirable as they are designed to han¬ 
dle strong local correlations, they are less accurate in 
treating non-local interactions, and are also limited in 
momentum-space resolution. Although it neglects strong 
correlation physics, the gRPA has the advantage that it 
treats local and non-local interactions on the same foot¬ 
ing, and is unbiased with respect to wavevector and to 
the unit cell-resolved current pattern. Within gRPA, the 
2-particle vertex function includes both exchange and di¬ 
rect interaction diagrams and hence also generates combi¬ 
nations of both (see Fig.[^, while Green functions remain 
unrenormalized. 

Following Ref.l^, we project the interactions onto a 
set of 19 basis functions ( 7 ),,^(k) in orbital and momen¬ 
tum space, leading to a 19 x 19 matrix equation for the 
effective interaction vertex r^^{q), where q = (q, w), and 
i,j label the basis functions. The basis functions and 
the interaction vertex are the same as in Ref.l^. Closing 
F*-^ (q) on the left and right with current vertex functions 
yields the susceptibility 


Xr'/ 5 / 3 '( 9 ) = xLw.( 9 ) -Y.A:^^^r{q)r^{q)A^fr{q), 

(16) 

where 



A 


i,Li 

OiOi' 


= “<q (k - q, q), 

k/iix 

= ^ E /iL'(k,q)M;l®'“<^(cc)5^,,(k) (17) 


Ai,Rl 

-^olol' 





,5y,(k)5,*,(k)^e>(k + q)^;^(k + q) (18) 

/(gk)-/K+q) 


I denotes bond indices, (k,q) are matrix elements 
of the current operators which are explicitly defined in 
Table [Tj 5'Qy(k) is the ath element of the vih eigenvector 
of Ho^), <5'*,,(k) is its complex conjugate, are the 
eigenvalues, f{E) is the Fermi function, and i and 5 in 
Eq. are the complex constant and a small broadening 
parameter respectively. The bare current susceptibility 
Xa ( 9 ) and the functions (q) and (q) differ 
from Ref.^ as they contain current operators. 


£ 

ee' 

^efl'(k,q) 

1 

dx 

^^^gh9x-l-fcx)/2 

1 

xd 


2 

dy 

_^^^gh9»+fcy)/2 

2 

yd 


3 

xy 


3 

yx 

,•+ „i(qx+k^-ky)/2 

Ol/ppC^ 

4 

dx 

/'2 

4 

xd 


5 

dy 


5 

yd 


6 

xy 

a p-i(qv+k:c-ky)/2 

6 

yx 

H f,-i(q!i:+kx-ky)/2 

LLppC, 

7 

xy 

nf p^(<?y+fcx+fcy)/2 

OLppC. 

7 

yx 

4i p—'('Jx+fex+fcy)/2 

tLppC. 

8 

xy 

j^^^g-*(ly-|-fex-|-fey)/2 

8 

yx 

nf pK^si+kx+ky)/2 

OLppC. 


TABLE I: Matrix elements of the current operator, i is the 
imaginary constant. The overall sign of each term results 
from three factors: the complex constant in the current oper¬ 
ator definitions, the sign of the hopping terms, and the gauge 
transformation. 


We search for the existence of spontaneous currents by 
following the evolution of the current susceptibility upon 
cooling. The instability is signalled by a divergence of 
the momentum-resolved susceptibility at zero frequency 


IV. RESULTS 

Figure shows typical results for the current sus¬ 
ceptibility. The inset shows that the matrix element 
xii{ci,uj = 0), corresponding to currents along the d-px 
bonds, becomes strongly peaked at q = Q = (tt, tt) as the 
temperature is lowered. This peak indeed diverges upon 
cooling to the critical temperature near T = 0.01 (main 
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FIG. 3: Diagrammatic structure of gRPA: (a) Interacting 
susceptibility, (b) vertex function, (c) effective interaction. 
Reprinted with permission frorrP^. ©American Physical So¬ 
ciety. 


panel), which signals an instability towards a current- 
carrying state. The ordering wavevector Q of the ttLC 
phase is the same as in the DDW scenario,*^ and should 
be contrasted with the On phase, ^ for which q = 0 . 

The example result in Fig.|^was obtained for Vpd = 2.2 
and Vpp = 1. However, the instability towards a ttLC 
phase persists when Vpp = 0, and hence is driven by the 
Coulomb repulsion Vpd between copper and oxygen or¬ 
bitals. In fact, also the local interactions Ud and Up have 
no effect on the ttLC instability. As we have explicitly 
verified, the staggered current instability originates from 
the exchange (ladder only) diagrams. 

As previously establishecP^, an ICO with a predomi¬ 
nant charge redistribution between Opx and Opy orbitals 
can be generated by the Coulomb repulsion Vpp between 
Opx and Opy orbitals. Based on the inter-orbital dis¬ 
tances, we expect Vpd > Vpp, which implies that loop 
currents emerge at higher temperatures than ICO. In¬ 
deed, for our parameter values, the critical temperature 
for the ttLC instability is about twice as large as the crit¬ 
ical temperature for ICO. 

To determine the bond-resolved ttLC pattern, we cal¬ 
culate the eigenvector of the leading eigenvalue of the 
current susceptibility matrix. In the current opera¬ 
tor basis Ji(q),..., J8(q), this normalized eigenvector 
is [0.48,-0.48,-0.15, 0.48,-0.48,-0.15,-0.15,-0.15]; 
this eigenvector reveals the direction and the relative 
magnitudes of the currents on the eight inequivalent 
bonds: all bonds are involved in the ttLC instability, and 
for the selected parameter set the currents along the p-d 
bonds are about three times stronger than those along 
the p-p bonds. The relative strength of the p-d and p-p 
currents varies with the ratio tpp/tpd- The wavevector Q 
of the instability further implies that the pattern alter¬ 
nates between adjacent unit cells. We thus obtain the 
cartoon shown in Fig. in which two distinct (green 
and black) and interpenetrating loop currents are evi¬ 
dent. This pattern is similar to the previously proposed 


xii(q,w = 0) 



FIG. 4: (Golor online) Temperature evolution of Xii(Q,t<J = 
0) for Q = (rr,7r). (inset) Xn(q,a; = 0) at T = 0.011. The 
parameter values are Vpd = 2.2, Vpp = 1, and the hole density 
is p = 0.13. 


current-carrying phases of either the DDW or the 
staggered flux-phase typd^Il*^, but differs in having two 
circulating current loops one of which involves oxygen 
orbitals only. 

The instability towards spontaneous ttLCs will natu¬ 
rally reconstruct the Fermi surface. To explore this we 
implement the ttLCs on the mean-field level. The starting 
point for this calculation is the general identity noUb = 
{J^hJnh/iih+ng+nb)/2 that is true for any pair of orbitals 
a and Hence, the non-local density-density interac¬ 
tions Can be decouplcd by intro¬ 

ducing the current amplitudes Ziajp = {Jiajis) = 

The mean-field version of the interorbital Coulomb inter¬ 
actions thus reads: 


^MF ~ + fly) + Cdfld 


E Via,jl3 f 


( 20 ) 

where ip = Vpd + 2V),p and id = 2Vpd renormalize the 
orbital energies. The intraorbital interactions lead to ad¬ 
ditional Hartree shifts of the orbital energies; these are 
assumed to be already included in and Cp. We obtain 
the mean field Hamiltonian Hmf = Ho + H'^^p, with 
Hmf = 'f'kHMF(k)4'k, = (^^k, 'I'k-t-q)'^, and 


HMF(k) 


Hi(k,Q) 


Ho(k) Hi(k,Q) 
Hl(k,Q) Ho(k + Q)J ’ 

0 %RpdS X 'iHpd^y 

RpdCx 0 ~fff‘ppl'x^y 
Rpdl-y ff'pp^x^^y ^ 


( 21 ) 

,( 22 ) 


where Rpd — ZpdVpd/tpd, Rpp — ‘^ZppVpp/tpp, — 
sm[{kx + Qx)/% c'x = cos[{kx + Qx)/2], and Sy,c'y are 
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FIG. 5: (Color online) Left: Normal (red) and reconstructed 
(blue) Fermi surfaces. The loop current strengths are set to 
Zpd = 0.04 and Zpp = ZpdjS. The hole density here is p = 0.10. 
Right: Density of states near the Fermi energy. 




defined accordingly, and the orbital energies in Hq (k) are 
shifted by id and ip. 

We show the reconstructed {zpd = 0.04, Zpp = Zpd/3) 
and normal {Zpd = Zpp = 0) Fermi surfaces and den¬ 
sities of states in Fig. For very small Zpd and Zpp, 
hole pockets, as well as small electron pockets near (tt, 0 ) 
and symmetry-related points, are formed. The shapes of 
the pockets depend on the strength of the loop currents, 
the hole filling, and the curvature of the unfolded Fermi 
surface. With increasing Zpd and Zpp, the electron pock¬ 
ets rapidly disappear, and only the hole pockets persist 
(Fig. g. If one takes tpd from band structure calcula¬ 
tions, Zpd = 0.04 corresponds to a p-d current of 16 /iA; 
empirical bandwidths are 3 times smallergiving a cur¬ 
rent of 5 pA and a plaquette magnetic moment in the 
range of O.OS/is to O.OO/is (see Appendix [ b|. 

We propose that the pseudogap appearing at T* 
arises from the ttLC shown in Fig. Physical 

propertie s of stagg ered currents have been discussed 
previousljEHillllSir and we focus here on aspects related 
to the recent discoveries of charge order and time-reversal 
symmetry breaking. The main distinguishing feature of 
the pseudogap is the depletion of spectral weight along 
regions of the Fermi surface near (±7r, 0) and (0, ±7r). 

This leads to a pseudogap in the density of states, as 
shown in Fig. for the ttLC phase. In our calculations 
Zpd is about one third of the peak-to-peak pseudogap; ex¬ 
perimental pseudogaps of ^ 100 me\® therefore suggest 
Zpd ~ 0.033eV, corresponding to a p-d current of ^ 7/xA, 
consistent with the estimate above. 


V. DISCUSSION 

Within the phase with staggered loop currents subse¬ 
quent phase transtions are likely to occur. Notably, we 
showed previously that persistent discrepancies between 
theory and experiment regarding the ordering wavevector 
q* of the charge ordered phase are resolved, if the charge 
order emerges from a preexisting pseudogap phase, rather 
than causing it. In Ref.^, a spin-density wave (SDW) 


with ordering wavevector Q = (tt, tt) was invoked ad hoc 
as a proxy for the pseudogap in underdoped cuprates. 
While the presence of a static SDW is not supported by 
experiment, we view the ttLC phase instead as a viable 
alternative phase out of which a charge-density wave will 
form with a q* that connects adjacent hole pockets. In 
the coexistence with charge order, the Fermi surface of 
the ttLC phase will further reconstruct. The Fermi sur¬ 
face in the coexistence phase should then serve as the ba¬ 
sis to analyze the quantum oscillation experiments which 
reported evidence for the existence of hole pockets.l^^lini 

The ttLC phase shares neutron-scattering signatures 
with the DDW state, specifically an elastic magnetic peak 
centered at Q = (tt, tt) and the opening of a spin excita¬ 
tion gap. Soon after the original proposal of the DDW 
state it was arguecP^^ that the neutron-scattering data 
obtained by Mook et al. in underdoped YBCO^ are 
consistent with the expected features of a DDW state. 
Other subsequent neutron scattering measurements on 
oxygen ordered ortho-II YBCO insteacP^ reported no 
evidence for the predicted characteristics of an ordered 
DDW state. The conflicting results of these experiments 
have remained unresolved. 

An obvious signature of the ttLC phase is that it breaks 
time-reversal symmetry. It does not, however, generate 
a polar Kerr effect because the pattern in Fig. pre¬ 
serves mirror symmetries. Given this result, an expla¬ 
nation for the observed nonzero Kerr angle at TKen- re¬ 
quires the onset of a further transition that eliminates 
these symmetries.^ This could naturally occur, for ex¬ 
ample, with the appearance of incommensurate charge 
order at T^o- The e xperimental doping dependences of 
Tkerr^ and 7).j 54 | 55 | however, different; a possible 
connection between the two is therefore not obvious. 

It is nevertheless possible that the different symmetries 
associated with the charge ordering transition are broken 
at distinct temperatures.122 Charge order involves both a 
continuous broken symmetry associated with the spatial 
lock-in of the charge modulation, and a discrete broken 
symmetry associated with the orientation of q* In this 
scenario, Txerr may signal the onset of a charge-nematic 
phase, in which only the discrete rotational symmetry is 
broken, while Tco may mark the pinning of the charge 
modulation by disorder.l^ 

We are led to trace the origin of the intriguingly com¬ 
plex physics of underdoped cuprates to the distinct phe¬ 
nomena which emerge from the Coulomb interactions in 
the Cu02 planes: the local Coulomb repulsion on the 
Cu d-orbital is the source of antiferromagnetism in the 
undoped compounds and of spin-fluctuation mediated d- 
wave superconductivity upon doping. From the results 
of this work we conclude that the non-local interaction 
Vpd can cause an orbital current instability, while the 
non-local interaction Vpp is responsible for the charge re¬ 
distribution between Opx and Opy orbitals and incom¬ 
mensurate charge order. Vpp and Vpd weaken the spin- 
fluctuation mediated pairing interaction, which suggests 
a possible reason why Tc < Tco below optimal doping. 
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The physics of underdoped cuprates therefore appears 
to reflect the mutual competition and/or coexistence of 
these ordering tendencies. 

VI. CONCLUDING REMARKS 

Within gRPA, we find no trace of a q = 0 current 
instability that has been proposed to exist in the three 
band model. Recently, Weber et a/. concluded that the 
charge-transfer energy, ej, — Cp, is one of the key param¬ 
eters for this instability. However, over an exhaustive 
range of this parameter, our calculations did not trace 
any instability or a sizable enhancement of the q = 0 
current susceptibility. Intra-unit cell magnetism has b een 
inferred from spin-flip neutron scattering experiments^. 
At this stage, our calculations offer no explanation for 
this observation. 

The close connection between axial charge order with 
a d-wave form factor!^ and oxygen orbitals makes the 
three-band model a necessary starting point for the mi¬ 
croscopic theory. In this context, the focus has so 
far been on spin fluctuations, and hence on the local 
Coulomb interactions Ud and Up. Our findings, however, 
point to the important role of non-local interactions. The 
development of accurate numerical tools capable of han¬ 
dling multi orbital models with intermediate to strong 
and non-local interactions is, therefore, on demand. 

We finally note that, for the typical parameter values 
used in this work, an SDW phase will, within gRPA, in 


fact set in at a higher temperature than the ttLC phase. 
This SDW is driven by the local Coulomb interaction Ud 
on the Cu sites, and it is known that correlation effects 
beyond the approximations discussed here suppress this 
SDW. Indeed, the feedback of spin fluctuations renor¬ 
malizes the electronic structure in a way that reduces 
the tendency towards antiferromagnetism. We have ne¬ 
glected this aspect here, since we are primarily interested 
in the charge degree of freedom in the underdoped part of 
the cuprate phase diagram where pseudogap and charge 
order occurs. Furthermore, we confirmed that the lo¬ 
cal interactions Ud and Up, the key parameters which 
control antiferromagnetism and spin fluctuations, do not 
have any effect on the emergence of the ttLC phase. 

It remains yet to be explored how self-energy correc¬ 
tions in our present approach will influence the instabil¬ 
ity towards spontaneous loop currents. This is on the 
agenda for further work on loop currents in the three- 
orbital model for cuprates. 


Acknowledgements 

A.P.K. and S.B. were supported by the DFG through 
TRR80. W.A.A. acknowledges support by the National 
Sciences and Engineering Research Council (NSERC) of 
Canada. We thank L. Chioncel for helpful discussions. 
The calculations were performed on the Linux Cluster of 
the LRZ in Carching. 


^ A. Shekhter, B. J. Ramshaw, R. Liang, W. N. Hardy, D. A. 
Bonn, F. F. Balakirev, R. D. McDonald, J. B. Betts, S. C. 
Riggs, and A. Migliori, Nature (London) 498, 75 (2013). 

^ B. Fauque, Y. Sidis, V. Hinkov, S. Paillies, C. T. Lin, 

X. Chaud, and P. Bourges, Phys. Rev. Lett. 96, 197001 
(2006). 

® Y. Li, V. Baledent, N. Barisic, Y. Cho, B. Fauque, Y. Sidis, 
G. Yu, X. Zhao, P. Bourges, and M. Greven, Nature 455, 
372 (2008). 

^ Y. Li, V. Baledent, N. Barisic, Y. C. Cho, Y. Sidis, G. Yu, 

X. Zhao, P. Bourges, and M. Greven, Phys. Rev. B 84, 
224508 (2011). 

® Y. Sidis and P. Bourges, J. Phys.: Conf. Ser. 449, 012012 
(2013). 

® C. M. Varma, Phys. Rev. B 73, 155113 (2006). 

^ C. M. Varma, J. Phys.: Condens. Matter 26, 505701 
(2014). 

® C. Weber, A. Lauchli, F. Mila, and T. Giamarchi, Phys. 

Rev. Lett. 102, 017005 (2009). 

® C. Weber, T. Giamarchi, and C. M. Varma, Phys. Rev. 
Lett. 112, 117001 (2014). 

R. Thomale and M. Greiter, Phys. Rev. B 77, 094511 
(2008). 

Y. F. Kung, C.-C. Chen, B. Moritz, S. Johnston, 
R. Thomale, and T. P. Devereaux, Phys. Rev. B 90, 224507 
(2014). 

J. Xia, E. Schemm, G. Deutscher, S. A. Kivelson, D. A. 


Bonn, W. N. Hardy, R. Liang, W. Siemons, G. Koster, 
M. M. Fejer, and A. Kapitulnik, Phys. Rev. Lett. 100, 
127002 (2008). 

R. -H. He, M. Hashimoto, H. Karapetyan, J. D. Koralek, 

J. P. Hinton, J. P. Testaud, V. Nathan, Y. Yoshida, H. Yao, 

K. Tanaka, W. Meevasana, R. G. Moore, D. H. Lu, S.- 
K. Mo, M. Ishikado, H. Eisaki, Z. Hussain, T. P. Dev¬ 
ereaux, S. A. Kivelson, J. Orenstein, A. Kapitulnik, and 

Z.-X. Shen, Science 331, 1579 (2011). 

Y. Wang, A. Chubukov, and R. Nandkishore, Phys. Rev. 
B 90, 205130 (2014). 

T. Wu, H. Mayaffre, S. Kramer, M. Horvatic, C. Berthier, 
W. N. Hardy, R. Liang, D. A. Bonn, and M.-H. Julien, 
Nature (London) 477, 191 (2011). 

G. Ghiringhelli, M. Le Tacon, M. Minola, S. Blanco- 
Canosa, C. Mazzoli, N. B. Brookes, G. M. De Luca, 
A. Frano, D. G. Hawthorn, F. He, T. Loew, M. M. Sala, 

D. C. Peets, M. Salluzzo, E. Schierle, R. Sutarto, G. A. 
Sawatzky, E. Weschke, B. Keimer, and L. Braicovich, Sci¬ 
ence 337, 821 (2012). 

J. Chang, E. Blackburn, A. T. Holmes, N. B. Christensen, 
J. Larsen, J. Mesot, R. Liang, D. A. Bonn, W. N. Hardy, 
A. Watenphul, M. v. Zimmermann, E. M. Forgan, and 

S. M. Hayden, Nature Phys. 8, 871 (2012). 

E. H. da Silva Neto, P. Aynajian, A. Frano, R. Comin, 
E. Schierle, E. Weschke, A. Gyenis, J. Wen, J. Schneeloch, 

Z. Xu, S. Ono, G. Gu, M. Le Tacon, and A. Yazdani, 



7 


Science 343, 393 (2014). 

R. Comin, A. Frano, M. M. Yee, Y. Yoshida, H. Eisaki, 
E. Schierle, E. Weschke, R. Sutarto, F. He, A. Soumya- 
narayanan, Y. He, M. Le Tacon, I. S. Elfimov, J. E. Hoff¬ 
man, G. A. Sawatzky, B. Keimer, and A. Damascelli, Sci¬ 
ence 343, 390 (2014). 

Y. Kohsaka, C. Taylor, K. Fujita, A. Schmidt, C. Lupien, 
T. Hanaguri, M. Azuma, M. Takano, H. Eisaki, H. Takagi, 

S. Uchida, and J. C. Davis, Science 315, 1380 (2007). 

W. D. Wise, M. C. Boyer, K. Chatterjee, T. Hondo, 

T. Takeuchi, H. Ikuta, Y. Wang, and E. W. Hudson, Nat. 
Phys. 4, 696 (2008). 

M. J. Lawler, K. Fujita, J. Lee, A. R. Schmidt, Y. Kohsaka, 

C. K. Kim, H. Eisaki, S. Uchida, J. C. Davis, J. P. Sethna, 
and E.-A. Kim, Nature (London) 466, 347 (2010). 

T. Wu, H. Mayaffre, S. Kramer, M. Horvatic, C. Berthier, 
W. N. Hardy, R. Liang, D. A. Bonn, and M.-H. Julien, 
Nat. Commun. 6 (2015). 

Y. Wang and A. Chubukov, Phys. Rev. B 90, 035149 
(2014). 

P. A. Lee, Phys. Rev. X 4, 031017 (2014). 

D. F. Agterberg, D. S. Melchert, and M. K. Kashyap, Phys. 
Rev. B 91, 054502 (2015). 

K. B. Efetov, H. Meier, and C. Pepin, Nat. Phys. 9, 442 
(2013). 

C. Pepin, V. S. de Carvalho, T. Kloss, and X. Montiel, 
Phys. Rev. B 90, 195207 (2014). 

D. Chowdhury and S. Sachdev, Phys. Rev. B 90, 134516 
(2014). 

L. E. Hayward, D. G. Hawthorn, R. G. Melko, and 
S. Sachdev, Science 343, 1336 (2014). 

A. J. Achkar, F. He, R. Sutarto, C. McMahon, 

M. Zwiebler, M. Hucker, G. D. Gu, R. Liang, D. A. 
Bonn, W. N. Hardy, J. Geek, and D. G. Hawthorn, 
arXiv: 1409.6787 (unpublished). 

R. Comin, R. Sutarto, F. He, E. H. da Silva Neto, L. Chau- 
viere, A. Frano, R. Liang, W. N. Hardy, D. A. Bonn, 
Y. Yoshida, H. Eisaki, A. J. Achkar, D. G. Hawthorn, 

B. Keimer, G. A. Sawatzky, and A. Damascelli, Nat. 
Mater. 14, 796 (2015). 

M. H. Fischer and E.-A. Kim, Phys. Rev. B 84, 144502 

( 2011 ). 

S. Bulut, W. A. Atkinson, and A. P. Kampf, Phys. Rev. B 
88 (2013). 

W. A. Atkinson, A. P. Kampf, and S. Bulut, New J. Phys. 
17, 013025 (2015). 

S. Chakravarty, R. B. Laughlin, D. K. Morr, and C. Nayak, 
Phys. Rev. B 63, 094503 (2001). 

S. Chakravarty, H.-Y. Kee, and K. Volker, Nature (Lon¬ 
don) 428, 53 (2004). 

V. J. Emery, Phys. Rev. Lett. 58, 2794 (1987). 

M. S. Hybertsen, M. Schliiter, and N. E. Christensen, 
Phys. Rev. B 39, 9028 (1989). 

P. B. Littlewood, C. M. Varma, S. Schmitt-Rink, and 

E. Abrahams, Phys. Rev. B 39, 12371 (1989). 

T. C. Hsu, J. B. Marston, and 1. Affleck, Phys. Rev. B 43, 
2866 (1991). 

D. A. Ivanov, P. A. Lee, and X.-G. Wen, Phys. Rev. Lett. 
84, 3958 (2000). 

K. Tsutsul, D. Poilblanc, and S. Capponi, Phys. Rev. B 
65, 020406 (2001). 

Y. He, J. Moore, and C. M. Varma, Phys. Rev. B 85, 
155106 (2012). 

K. Pasanai and W. A. Atkinson, Phys. Rev. B 81, 134501 


( 2010 ). 

R. B. Laughlin, Phys. Rev. B 89, 035134 (2014). 

I. Dimov, P. Goswami, X. Jia, and S. Chakravarty, Phys. 
Rev. B 78, 134529 (2008). 

T. Kurosawa, T. Yoneyama, Y. Takano, M. Hagiwara, 

R. Inoue, N. Hagiwara, K. Kurusu, K. Takeyama, 

N. Momono, M. Oda, and M. Ido, Phys. Rev. B 81, 094519 

( 2010 ). 

N. Doiron-Leyraud, C. Proust, D. LeBoeuf, J. Levallois, 

J. -B. Bonnemaison, R. Liang, D. A. Bonn, W. N. Hardy, 
and L. Taillefer, Nature (London) 447, 565 (2007). 

S. E. Sebastian, N. Harrison, and G. Lonzarich, Rep. Prog. 
Phys. 75, 102501 (2012). 

S. Chakravarty, H.-Y. Kee, and C. Nayak, Int. J. Mod. 
Phys. B 15, 2901 (2001). 

H. A. Mook, P. Dai, and F. Dogan, Phys. Rev. B 64, 
012502 (2001). 

C. Stock, W. J. L. Buyers, Z. Tun, R. Liang, D. Peets, 

D. Bonn, W. N. Hardy, and L. Taillefer, Phys. Rev. B 66, 
024505 (2002). 

M. Hucker, N. B. Christensen, A. T. Holmes, E. Black¬ 
burn, E. M. Forgan, R. Liang, D. A. Bonn, W. N. Hardy, 

O. Gutowski, M. V. Zimmermann, S. M. Hayden, and 
J. Chang, Phys. Rev. B 90, 054514 (2014). 

S. Blanco-Canosa, A. Frano, E. Schierle, J. Porras, 

T. Loew, M. Minola, M. Bluschke, E. Weschke, B. Keimer, 
and M. Le Tacon, Phys. Rev. B 90, 054513 (2014). 

®® L. Nie, G. Tarjus, and S. A. Kivelson, Proc. Nat. Acad. 
Sci. Ill, 7980 (2014). 


Appendix A: Current operators 

The current operator is conventionally defined as 

Jij =-iUj{clc^ - clci). (Al) 

At this point, however, it is important to take into ac¬ 
count that we are working with a Hamiltonian which is 
gauge transformed according to Equations In the 

transformed operator basis, the current operators in real 
space take the following form: 


Jid^jx 

= + c]a;Cid) 

(A2) 

Jix,jd 

= —tix,jd{c\iiCjx + C^jx^id) 

(A3) 

Jid.jy 

= '^id,jy{^id^jy + CjyCid) 

(A4) 

Jiyjd 


(AS) 

Jix,jy 

= ~Hix,jy{^ix^jy ~ ^jyCix) 

(A6) 

Jiy,jx 

= ~'itiy,jx{c\yCjx — C^jx^iy)- 

(A7) 


Appendix B: Estimating the resulting magnetic 
moment 

Here, we perform a simple estimate of the magnetic 
moment associated with the ttLC loop currents. 

We start by expressing the currents in absolute 
units. The unit of the current operators is = 


C{eV)l{eV ■ s) = Cjs. For z^d = 0.04ipd, the cur¬ 
rent along 'p-d bonds is l^d = x (e/?i) = 0.04tpd x 
(2.4 X = 9.6tpdpA, where tpd is measured in 

eV. Band structure calculations suggest tpd « 1.6 eV, 
yielding Ipd = 16 /iA; experimental bandwidths, how¬ 
ever, are typically a factor of 3 smaller than predicted by 
band structure calculations suggesting Ipd 5 /iA. This 
is comparable to other estimates of loop current ampli¬ 
tudes: within a cluster calculation of multi-orbital t-J 
model, the upper bound of 0//-like loop currents was pre¬ 
viously estimated to be between 5pA and 15/iA for differ¬ 
ent parameter seti^^ similarly, in the single band DDW 
studies, the staggered loop currents were estimated to be 
~ 7pA by assuming a DDW gap value of ^ 0.03e\^. 

Next, we calculate the magnetic moments of loop cur¬ 
rents using the formula M = Irj where I is the current, 
and r] is the area enclosed by the loop. This is indeed 
a crude calculation, however we believe this should yield 
a qualitatively correct number. The total magnetic mo¬ 


ment in a given plaquette has two contrubutions: 
and IgTjg due to two independent current loops (shown 
in black and green in Fig. 1). We set /(, = 5/rA, and 
since it circulates around the whole plaquette, r]b = a^. 
As given by the eigenvectors of Ig ~ Ipd /and 

it circulates around an area of pg = a^/2. Thus, the 
magnetic moment is calculated as 


M = IbVb + IgVg (Bl) 

= Ipd X + Ipd/S X a^/2 (B2) 

= 7VaV6 (B3) 

= 7(5/rA)(3.85 x 10 "^°to)V6 (B4) 

« OmpB- (B5) 


However, since Ipd is shared by the neighbouring plaque- 
ttes, the effective moment for individual plaquettes might 
be reduced to ~ O.ObpB- 


